Electron Transport Driven by Nonequilibrium Magnetic Textures 
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Spin-polarized electron transport driven by inhomogeneous magnetic dynamics is discussed in the 
limit of a large exchange coupling. Electron spins rigidly following the time-dependent magnetic 
profile experience spin-dependent fictitious electric and magnetic fields. We show that the electric 
field acquires important corrections due to spin dephasing, when one relaxes the spin-projection 
approximation. Furthermore, spin-fiip scattering between the spin bands needs to be taken into 
account in order to calculate voltages and spin accumulations induced by the magnetic dynamics. 
A phenomenological approach based on the Onsager reciprocity principle is developed, which allows 
us to capture the effect of spin dephasing and make a connection to the well studied problem of 
current-driven magnetic dynamics. A number of results that recently appeared in the literature are 
related and generalized. 
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Interest in magnetic heterostructures/ which was ini- 
tially fueled by the discovery of the giant magnetoresis- 
tance and, a decade later, by the current-induced switch- 
ing in spin valves and related systems, has more recently 
spilled over into current-driven phenomena in magnetic 
bulk, individual magnetic films, and nanowires.~ Par- 
ticular attention was given to the problems of current- 
driven Doppler shift of spin waves, magnetic instabilities, 
and domain-wall motion. The latter has also enjoyed a 
very vibrant experimental activity, which is in part mo- 
tivated by a promising application potential in spintron- 
ics. The past yeari^i^i^ saw a revival of interest in the 
inverse effect of electromotive forces induced by the time- 
dependent magnetization, which were previously studied 
in various physical contexts (see, e.g., Refs.SIllH)- this 
paper, we will exploit the reciprocal relation between the 
two phenomena, which will allow us to understand im- 
portant spin-dephasing corrections to the electromotive 
force. Such corrections were first mentioned in Ref. and 
the Onsager principle in the present context was invoked 
in Ref. H. Reference H reported the magnetically- induced 
electromotive forces as a manifestation of the position- 
dependent Berry phase accumulation and Ref. J) consid- 
ered these forces acting on semiclassical wave packet mo- 
tion, mainly reproducing results from Ref. For com- 
pleteness, it should also be mentioned that a much earlier 
paperlS already contains some seminal phenomenological 
insights related to the problem of the electric response to 
the magnetic domain- wall dynamics. 

In the following, we start by recalling how most of the 
results recently discussed in the literature can be cap- 
tured by an SU(2) gauge transformation together with 
the projection of spins on the magnetic directioni^ The 
corrections due to the remaining transverse spin dynam- 
ics are governed by spin dephasing, which have already 
been studied for the reciprocal process of current-driven 
magnetic dynamics,^ and can be translated to the cur- 
rent problem by the Onsager principle. We will develop a 
general framework, that will allow us to relate and gener- 
alize the more specialized cases, which were recently stud- 



ied using different methods i^i^i^i^ Finally, we will derive 
spin-charge diffusion equations, accounting for spin-flip 
scattering and respecting local charge neutrality, which 
is necessary in order to relate the microscopic electromo- 
tive forces to measurable quantities, such as an induced 
voltage and spin accumulation. 

Most of our analysis will pertain to the following time- 
dependent Hamiltonian: 
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Here, is the contribution due to spin-relaxation pro- 
cesses, which will be characterized by a Bloch-type Ti 
spin flipping and T2 spin dephasing, and Vc(j,t) stands 
for a Hartree charging potential, which will be taken into 
account only insofar as enforcing local charge neutral- 
ity. Axe is the ferromagnetic exchange band splitting, & 
is the vector of Pauli spin matrices, and m stands for 
the local magnetization direction unit vector, so that the 
magnetization is given by M = Mm. The exchange field 
Axcin may in practice be provided by localized magnetic 
d orbitals (as in the so-called s — d model) or it may 
self-consistently be governed by the itinerant electron 
spin density (as in the Stoner model or local spin-density 
approximation) ^ We will first perform a microscopic cal- 
culation for the idealized Hamiltonian ^ neglecting Her 
and subsequently utilize the Onsager theorem to capture 
the spin-dephasing corrections. The spin- flip scattering 
will be included phenomenologically in the final diffusion 
equation. 

By disregarding Ho-, we can perform an SU(2) gauge 
transformation by rotating m to point along the z axis 
for all r and This is conveniently achieved by the 

Hermitian spin-rotation matrix U ^ & n (such that U = 
= U~^), where n is the unit vector n (x m + z. It is 
easy to see that U{a- ■ m.)U = az (since U corresponds 
to a TT-angle spin rotation around n). By applying this 
gauge transformation to the spinor wave function, we get 
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for the transformed Hamiltonian 

where the SU(2) vector potential is given by Ai = 
ihUViU — —ha ■ (n x Vin) and the SU(2) ordinary po- 
tential \sV — —ihUdtU — ha- ■ (n x dtn) (setting the par- 
ticle charge and speed of light to unity), p = —ifiV is the 
canonical momentum. If the exchange field Axe is large 
and the magnetic texture is sufficiently smooth and slow, 
we can project the fictitious potentials on the z axis as 
V V&z, where V — hz-{nxdtn) = hs\T?{9 /2)dt4'^ and 
similarly for the vector potential, Ai — ~hz ■ (n x Win) — 
—hsin^{9/2)\7i(j), where {0,<j)) are the spherical angles 
parametrizing m. We thus get for the effective electric 
fieldi^^^ 

E = -dtA - VI/ = ^ sin amiVcl)) - {M){V0)] , 

(3) 

or, written in the explicitly spin-rotationally-invariant 
form, 

E, = • {dtm X V,m) . (4) 

The effective magnetic field 

B = V X A = - sme{V(t)) X (V0) 

= -e'^^^iJiiiVrnk x Vnij) , 

where e^^^ is the antisymmetric Levi-Civita tensor, and 
a summation over repeated indices is implied. The total 
force on a spin-^ (J,) electron moving with velocity v is 
thus given by 

FT(i) = ±(E + v x B) . (6) 

Eqs. ^ and ([5]) were first derived in Ref.dand recently 
rederived within a semiclassical wave-packet analysis.— 
The gauge-transformation based approach^ puts the re- 
sult into a broader perspective, allowing us, for exam- 
ple, to consider the effect of the magnetic field (O on 
the quantum transport corrections, such as a weak local- 
ization, as well as to include spin-independent electron- 
electron interactions, which would not modify fictitious 
fields ([3]) and (|5]). Note, in particular, that the magnetic 
field ([5]) can in practice be quite large: For example, for 
a static magnetic variation on the scale of 10 nm, the 
corresponding fictitious field is of the order of 10 T. It 
is most convenient to estimate the strength of the elec- 
tric field ([21) by the characteristic voltage it induces over 
a region where the magnetization direction fiips its di- 
rection: hw / where u> is the frequency of the magnetic 
dynamics. In the following, we will concentrate on the 
semiclassical spin and charge diffusion generated by the 
effective electric field E. In order to make a closer con- 
nection to the experimentally relevant quantities, we will 



need to take into account spin relaxation and also enforce 
local charge neutrality for electron diffusion. 

The role of spin relaxation can be twofold. First of all, 
spin accumulation, which will generally be generated by 
the spin-dependent force ([5]) will relax, characterized by 
the longitudinal spin- flip time Ti. There is also another 
more subtle effect, which is due to the dephasing of elec- 
tron spins following a dynamic magnetic profile, since the 
exchange field Axe is not infinite and spins do not per- 
fectly align with the local magnetization. Hence, there is 
generally a finite spin misalignment, which dephases with 
a characteristic time T2. This gives corrections to the re- 
sults obtained by a rigid projection of spins on the local 
magnetization direction. We will see that such correc- 
tions turn out to be important for the currents generated 
by magnetic dynamics, in the same sense that analogous 
corrections are crucial for understanding current-induced 
magnetic motioni^ 

Let us now take a step aside, by recalling the general 
expression for the dynamics of an isotropic ferromagnet 
well below the Curie temperature)^ 

dM = - 7M X Heff + — M X dtM 

,i<£L_£il,C)(i-^M.)v.M, 

which is valid for spatially smooth magnetic profiles (the 
so-called adiabatic approximation) and weak currents. 
Here, a is the Gilbert damping constant, /3 is another di- 
mensionless phenomenological parameter whose physical 
meaning will be discussed later, /i is the electrochemical 
potential, as is the spin-s conductivity (along the local 
magnetization direction m) relating particle currents to 
V^, S is the equilibrium spin density of the ferromag- 
net along m, and 7 = M/S is the gyromagnetic ratio. 
Recall that the effective field Hoff is the quantity defined 
to be thermodynamically conjugate to the magnetization: 
Hefi = (note the sign difference from the standard 

definition), where J- is the free energy and i9m stands 
for the functional derivative. The other thermodynamic 
variable we will consider is the electron density p{r,t), 
whose thermodynamically conjugate counterpart is the 
electrochemical potential fi = dpJ-. 

Suppose we perturb the electron density with respect 
to an equilibrium with some static magnetic texture 
and uniform chemical potential, and consider the en- 
suing magnetic response. Eq. ([7]) then describes the 
nonequilibrium coupling of the magnetization dynam- 
ics to the electron density's thermodynamic conjugate, 
which is slightly out of equilibrium. The Onsager reci- 
procity principlc^"^ allows us to immediately write down 
the response of the electron density to a small modulation 
of the effective field Hoff, with respect to an equilibrium 
configuration. To simplify things, let us for a moment 
disregard Gilbert damping a in Eq. ([7]) and return to 
include it later on. An electric response to a magnetic 
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perturbation then becomesi^ 
aip = -2:M^L^V,{Heff •[(l + /3mx)V,m]} . (8) 

By comparing Eq. ((S]) with the continuity equation dtp — 
—'^iji, we can identify the particle current as 

lh{a^^ "^^ Heff ■[(l + /3mx)V,m] . (9) 
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Since for each spin species, js = CTsFg, where is the 
effective force, we finally get for the latter Fj^j = ±F— , 
where 



= - (m X dtm) ■ [(1 + /3mx) V^m] 

= ^ [m • {dtm X Vim) + f3 {dtm ■ Vim)] 



(10) 



after inverting the magnetic equation of motion fT]) in 
order to express the effective field Hoff in terms of the 
magnetization dynamics m(r,t). (Note that since the 
currents themselves are now generated by the magneti- 
zation dynamics, we can neglect their backaction on the 
magnetic response, when inverting the equation of mo- 
tion to express Hoff in terms of m, since it would give 
rise to higher-order terms.) Equation (jlOp is a key result 
of this paper. It is also easy to show that taking into 
account Gilbert damping a has no consequences for the 
final result pO)) [after rewriting Eq. ([7]) in the Landau- 
Lifshitz form, in order to eliminate the dt term on the 
right-hand side and thus make the equation suitable for 
the Onsager theorem]. This is not surprising, since the 
physics of the Gilbert damping a does not have to be 
related to the magnetization — particle-density coupling 
that determines the force PTI)) .^ 

Physically, the /3 correction in Eq. PH]) is related 
to a slight spin misalignment of electrons propagating 
through an inhomogeneous magnetic texture with the lo- 
cal direction of the magnetization m. In the limit of 
Axe — *■ oo, this misalignment vanishes and so should /3, 
reducing the result (fTU|) to Eq. (H]). Indeed, a microscopic 
derivation of Eq. ([7]) shows /3 ^ ?i/(T2Axc), where T2 is 
the characteristic transverse spin relaxation timci^ The 
P term in Eq. ((TO)) can thus be viewed as a correction to 
the topological structure of the electron transport rigidly 
projected on the magnetic texture, due to the remaining 
transverse spin dynamics and dephasing. Such a /3 cor- 
rection was first reported in Ref. i^, which used a very 
different and more technical language and did not bene- 
fit from the reciprocity relation with the current-driven 
magnetic dynamics ((T]). Our phenomenological deriva- 
tion of Eq. (fTO|) based on the Onsager theorem provides 
a much simpler framework for studying these subtle spin- 
dephasing effects. 

Let us now discuss the measurable consequences of 
Eq. ((T0|) in two simple scenarios sketched in Fig. [TJ Con- 
sider a nontrivial one-dimensional magnetic profile along 
the X axis, such as a magnetic domain wall in a nar- 
row wire, with negligible transverse spin inhomogeneities. 
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FIG. 1: Two simple scenarios for voltage generation by the 
magnetic dynamics: (1) Magnetic texture m{x,t), such as a 
domain wall along the x axis, is steadily rotating around the 
X axis and (2) the same texture rigidly sliding along the x 
axis. In the former case, the force Fx acting on electrons is 
proportional to the frequency of rotation w, with the domi- 
nant term having a purely geometric meaning in terms of the 
position-dependent Berry-phase accumulation rate. (An al- 
ternative physical explanation can also be provided by trans- 
forming to the rotating frame of reference and applying the 
Larmor theorem.) In the case of the sliding dynamics, the 
leading contribution to the magnetically-induced force is pro- 
portional to the spin-dephasing rate (parametrized by /?) and 
the "curvature" of the texture profile {dxtn}^. 



First, let us look at a steady rotation of the entire one- 
dimensional texture around the x axis, with a constant 
frequency uj. Then, dtm = ujx x m and 



dxFx = — —X ■ / (dm -I- ^m x dm) . (11) 



In the absence of spin dephasing f3, this result can be eas- 
ily understood by transforming into the rotating frame 
of reference: By the Larmor theorem, this corresponds 
to a fictitious field along the x axis: H' = —{fiuj/2)ax- 
For spins up (down) projected on the local magnetiza- 
tion direction, this corresponds to the potential V — 
^{fiu!/2)x ■ m. It is equally straightforward to interpret 
this result in terms of the rate of the Berry phase accu- 
mulation by spins adiabatically following the steady ex- 
change field precession;2iii which is proportional to the 
position-dependent solid angle enclosed by spin preces- 
sion. The (3 term in Eq. pT|) gives a correction to these 
idealistic considerations, which depends on the geometry 
of the magnetic texture. Next, we consider the voltage 
induced by a rigid translation of a one-dimensional mag- 
netic texture m{x—vt) along the x direction with velocity 
V. The corresponding force 



F,^~^f3v {d,mf 



(12) 



is then entirely determined by the /3 term, which drags 
spins down along the direction of the magnetic texture 
motion and spins up in the opposite direction. This is 
analogous to the current-driven domain wall velocity in 
one dimension, which, for smooth walls and low currents, 
is proportional to 
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Finally, we need to include spin-flip relaxation time 
Ti and derive spin-charge diffusion equations, enforcing 
local charge neutrality. Assuming diffusive transport, the 
force pUj) can now be added as a contribution to the 
gradient of the effective electrochemical potential. The 
diffusion equation for spin-s particles is then given by 

{dt - DsW^)ps +as{sV-F- V Vc) = ^ - ^ , (13) 

where ps is the nonequilibrium (spin-s) particle density, 
Ds is the diffusion coefficient, and is the spin-flip time. 
Recall that the conductivity is related to the density of 
states Ns by the Einstein's relation: as = NgDs- Vc is the 
electric potential, which has to be found self-consistently 
by enforcing local charge neutrality. Note that the equi- 
librium considerations require that Ts/T_g = Ns/N^^. 
We should also stress that the force PH)) may have a finite 
curl, so that we cannot generally describe it by a ficti- 
tious potential. After straightforward manipulations, we 
can decouple the diffusion equation for the spin accumu- 
lation fia (defined as the difference between the spin-up 
and spin-down electrochemical potentials, divided by 2) 
from the average electrochemical potential /i as follows: 

(14) 

-V^/i = P (V^/i, - V • F) . 
Here, P = ((T| — crx)/(CT| + cr^) is the conductivity polar- 



ization, D = {D^+Di)/2~P{D^-Di)/2 is the effective 
spin-diffusion constant, and = t^^ ^'''"i"^ char- 
acteristic T^^ rate for spin flipping. If V x F = 0, we 
can integrate the second equation to express the electro- 
chemical potential gradient in terms of the force F and 
the spin accumulation gradient as follows: 

= P(F - Vm.) , (15) 

assuming the appropriate boundary conditions. Accord- 
ing to Eqs. (|14p . the spin accumulation decays in the 
absence of the force F on the scale of the spin-diffusion 
length Asd = \J Dt . Away from the dynamic magnetic 
texture (on the scale of Asd), the generated electrochemi- 
cal potential (|15p will then be determined simply by inte- 
grating the force F. In general, however, especially when 
V X F 7^ 0, one has to revert to Eqs. p^ . 

In summary, we theoretically studied electron trans- 
port generated by a dynamic magnetization texture. We 
reproduced and generalized the results that recently ap- 
peared in literaturC)^!^ revealing an intricate connec- 
tion with the theory of the current-induced magnetiza- 
tion dynamics 1^ We expect that in practice it is consider- 
ably simpler to solve this reciprocal problem, especially 
for including subtle corrections to the topological Berry 
phase structure of spins assumed to rigidly follow the 
time-dependent magnetic profile. 
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Note that the Onsager reciprocity relations are somewhat 
special in the present case; since one of the thermody- 



namic quantities (namely, the magnetization) is odd un- 
der time reversal, the matrix that couples the rates of the 
relaxation of the involved quantities to their thermody- 
namic conjugates is antisymmetric, which is subject also 
to changing the sign of the magnetic field and the equi- 
librium magnetizatioufi"^ Hence the overall minus sign in 
Eq. (|8]) and the sign change in front of 13. 
Strictly speaking, the analysis based on Eq. (O can only 
capture the total particle current ® . In order to rigorously 
calculate the spin-resolved forces, we would have to explic- 
itly take into account one more thermodynamic variable, 
namely, the nonequilibrium spin accumulation (which is 
proportional to the chemical potential mismatch between 
the up and down electrons along the local magnetization) 
and consider its action on the magnetic dynamics. The 
latter program would push us too much off track, and we 
choose not to pursue it here. We only wish to note that 
the P term describing the spin-dephasing correction to the 
fictitious field (|4]) may in general become different for the 
two spin species when the ferromagnetic exchange energy 
is comparable to the Fermi energy. 



